%% The following script is used to get F107 data from
% swpc website

icount =  1; 
for i = datenum(2011,8,1):datenum(2011,10,13),
s = urlread(['http://www.swpc.noaa.gov/ftpdir/warehouse/2011/RSGA/' datestr(i,'yyyymmdd') 'RSGA.txt']);
k = findstr(s,'IV.  Penticton 10.7 cm Flux');
f107(icount)=str2num(s(k+53:k+56));icount = icount + 1;
end;